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Abstract 

We have performed a partial-wave analysis of /("^-nucieon scattering in the momentum range from 
to 1.5 GeV/c addressing the uncertainties of the results and comparing them with several previous 
analyses. It is found that the treatment of the reaction threshold behavior is particularly important. 
We find a 2^=0 scattering length which is not consistent with zero, as has been claimed by other 
analyses. The T—Q phase shifts for £ > are consistent with a pure spin-orbit potential. Some 
indications for the production of a T—Q pentaquark with spin-parity D5/2+ are discussed. 



I. INTRODUCTION 

The interaction of the meson with the nucieon has a number of interesting features. Because 
of lack of annihilation of the meson anti-quark there is no three-quark intermediate state possible so 
3-q resonances are not possible. This lack of 3-body states leads to a feeble interaction, among the 
smallest among the strong forces. 

This fact has been useful forprobing the possible changes in nucieon structure and/or K+A^ inter- 
action in the nuclear medium [l|, [3, S fi 0, S- Following the measurement of the ratio of total 
cross sections [3i lEli a number of suggestions were made to explain the results, among the principal 
ones being a partial deconfinement of quarks, an interaction with exchanged mesons in the nucleus 
and a modification of the Tf+iV interaction through exchanged mesons. In order to understand which 
of these possible scenarios might be the right one, a detailed comparison with experiment is in order. 
One element in making calculations of the multiple scattering corrections which are crucial for this 
comparison is the availability of reliable phase shifts. 

Of equal importance is the use of the TfA^ system to test fundamental theories of hadronic interaction. 
There have been a number of theoretical studies of the N system using a variety of approaches 
d, [13, [O, m, [H, The phase shifts from K'^N scattering provide the principal data with which 
they compare. 

There have been very little data taken (but see Ij,j5|) in an ordinary sense since the most recent 
partial- wave analyses [ll, [11, [13, [11, [11 S [III Hi [p, [13, [H, [H, [ij [2|. The suggestion of the 
existence of a possible pentaquark resonance (see Ref. [2!l] for a list of experiments claiming to see it, 
some of which have been withdrawn) might change the point of view of any analysis. The existence 
of the pentaquark is now considered doubtful by many. 

A partial- wave analysis is particularly interesting in the region of the reaction threshold. The 
production of pions in the collision of two nucleons is often calculated with the "re-scattering" diagram 
which describes the formation of the A^'^ resonance by the exchange of one pion. In the case of 
scattering this simple re-scattering process does not exist since the cannot emit a single pion, and 
other mechanisms must be considered. The excitation by exchange of a p meson might be a natural 
process to consider, at least from the point of view of the exchange of heavy mesons. Additional 
information on this process can be obtained from a phase shift analysis which indicates which partial 
waves are participating in the production. 

Section II describes our fitting method. Section III gives a description of the T=l analysis. Section 
IV discusses the results of the T=l analysis, and Section V gives the method for the T=0 analysis. 
Section VI discusses the results of the T=0 analysis and Section VII gives an overall summary and 
discussion of the results. 
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II. FITTING PROCEDURE 



In the heart of any ampUtude analysis Ues the search for the (or a) minimum in a (or similar) 
measure for the best fit to the data. It is normally assumed that the numerical procedure for this 
minimization is straightforward and does not pose any problems, but that may not be the case. The 
technique that wc use is very pedestrian but seemingly very sure. The is minimized on each 
parameter in turn sweeping through a significant number of them (88 for the T=l case and 61 for 
the r=0 case including data normalizations). To minimize x^, each parameter is stepped by a fixed 
interval until the value of x^ increases. At this point in the search three values of x^ arc known at 
three values of the parameter. A parabola is then passed through the three points and the position of 
the minimum is predicted. This procedure, in some form, is common to most methods of minimization 
although several methods treat the full parameter space as a vector. There is a difficulty that arises 
due to the fact that x^ a-s a function of the parameter is very often not a perfect parabola. This 
means that the prediction of the minimum position (and hence x^ value) is not the true minimum in 
this region and the value of x^ predicted may exceed the one at the central (lowest) point. Since the 
middle point is very often the value obtained in the search on the previous parameter in the sequence, 
if the predicted value is always accepted, the "current minimum" x^ will increase in some cases. In 
our algorithm we test the predicted value of x^ against the central value and if it is greater, the central 
(previous best) value is used instead. 

This incorrect prediction is not a rare occurrence. Wc observed that it happened about 4% of the 
time when the search was far from the minimum and up to 40%-50% of the time when it was close 
to the minimum. If one decreases the step size, the deviation from a parabolic shape can be lessened, 
but there is a limit to this process since the difference between values of x^ needed to predict the new 
minimum becomes small compared with their values and with finite precision another source of error 
becomes important, even with double precision (which we use). Hence, there is an inherent limit to 
how well the minimum can be found. 

This limit depends on the details of the method used, of course, but also on the parameterization 
of the phase shifts. Different paramcterizations will have different functions to replace the parabola, 
or perhaps more practically, will have a different importance of third order terms. 

This inability to find perfectly the minimum (or minima) translates to a dependence of the final 
result on the starting values. Tests of the dependence on starting point were made by perturbing 
some of the parameters (typically two or three) until the x^ was very large. For x^ of the order of 
100,000 the search usually was unable to find a sensible minimum but for x^ in the range 10,000 to 
20,000 a minimum was found near, but not identical, to the principal fit. A large number of sweeps 
through the parameters was needed (several tens of thousands). The results of these tests are given 
in the discussions of the T=l and T=0 analyses. 

A x^ corresponding to each normalization was included in the total x^ of the fit by adding 



for each normalizing parameter, N, to the x^ coming from the individual data points. The value of 
0.03 was chosen for A7V since it is a typical value for normalization errors. 

The philosophy behind this procedure is that the normalization should be treated as an independent 
data point. In a typical model experiment the number of counts is registered in a set of counters and 
then those counts are multiplied by a normalization factor determined by the effective number of 
target particles, the beam flux etc. The errors in this normalization factor will be common to the 
cross sections for every counter. While this would appear to be the right statistical procedure, some 
care is needed in the interpretation of the results. Even if a significant discrepancy occurs in the 
normalization, it can be hidden if one simply looks at the total per data point in the case that 
there are a large number of points which accompany a single normalization. The large number of 
points dilute the total x^ so that it is only slightly larger than the number of degrees of freedom. A 
distinction based solely on total requires the separation of a value of the reduced x^ closer and 




(1) 
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closer to unity from unity. Practical difficulties arise since one can assume that the experimental 
errors have been estimated correctly only up to a certain level (experience would suggest something of 
the order of 5%). Hence the square of the errors in the denominator of the definition of leads to a 
sufficient error in its value that this method is impractical. This does not mean that a problem with the 
normalizations cannot be found, however. It simply means that one must examine the normalizations 
separately to verify that the the from the normalizations alone is not excessive. 



III. K+ PROTON (r=l) ANALYSIS 



The T=l data base used consists of 1880 data points including 1501 angular distribution points, 91 
total cross section points, 265 polarization points and 23 reaction cross section points. Much of these 
data are available from compilations 3 ^ 31 1 33 | but we try to give references to the original data as 

much as possible [H [H, [H, [H, [13, [aaJlMiil, US, S, H E m, [13, m, [ii, S [m, y, m . 

Of these points, four were removed because they were outliers. These were: two points from Ref. 

si at 0.686 and 0.717 GeV/c, and two from Ref. [H] at 0.713 and 1.029 GeV/c. The error on Ref. 

33| polarization at 1.330 GeV/c and cosO = 0.242 was doubled. The errors were doubled on the 
polarization set at 1.430 GeV/c [s^l since the variation from point to point is greater than the errors 
quoted. 
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FIG. 1: Distribution of the normalization factors for isospin 1. 



Partial waves through F7/2 were considered. The G and _ff waves were not needed to get a good 
fit. We write the partial- wave amplitudes as 



(S£±-l)e^^ 2i5i,± fcy\ 

i't± — , 0£± — rii±e '± {I) 

where the sign ± corresponds to j = £ ± i and ai is the Coulomb phase shift. No "inner" Coulomb 
corrections were made, the Coulomb effect being very small for the energies considered here. 
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LJ 


gfl (GeV/c) 


71 


72 


73 


Equation 


Sl/2 


0.4134 


0.0858 


0.5030 




El 


Pl/2 


0.3858 


0.1811 


-0.0163 




m 


P3/2 


0.4809 


1.0794 


2.1643 


-3.9024 


m 


D3/2 


0.6758 


2.1264 






m 


D5/2 


0.4304 


0.3604 


-0.2723 




m 


F5/2 


0.5286 


0.1778 






m 



TABLE I: Parameters for the representation of the inelasticity, r;^j for T=l using the form of Eq. [7] except 
for the Z)3/2 wave where Eq. [8] is used. 



LJ 


a (GeV/c)-(^^+i) 


6i,tt (GeV/c)-' 


&2,tt (GeV/c)-^ 


fos.tt (GeV/c)-'' 


Equation 


Sl/2 


-1.562 


-1.108 


0.217 




HI 


Pl/2 


-12.002 


31.139 






m 


P3/2 


13.357 


126.676 


-666.951 


1276.123 


n 


LJ 


c± {GeV/c)-'' 


d± (GeV/c)"^ 


e± (GeV/c)"^ 




Equation 


D3/2 


-2.984 


4.119 






El 


D5/2 


-1.702 


6.571 


-7.462 




m 


LJ 


/± (GeV/c)-^ 








Equation 


F5/2 


-0.415 








El 


FT/2 


0.089 








m 



TABLE II: Parameters for the representation of the phase shifts for T=l using the form of Eq. |4]for the S 
and P waves, Eq. [5] for the D waves and Eq. [6] for the F waves. 



The differential cross section, polarization and total cross section are expressed in terms of the 
amplitudes in the standard manner (see e.g. Hashimoto [l5| ) and will not be repeated here. The 
reaction cross section (essentially pion production) can be written as 

afl = ^ - r,j_) + (£ + 1)(1 - r,j^)] (3) 

SO is independent of the phase shifts and can be particularly important in determining the rn±. Here, 
and in what follows, q denotes the momentum in the center of mass (in GeV/c unless otherwise noted) 
and qf is q/hc, i.e., in femtometers. 

We have taken for the S- and P-wave phase shifts the form 



5i± = tan ^ 



1 + b^,i±q'' 



(4) 



The index, j runs from 1 to 2 for the S wave and 1 to 3 for the P3/2 wave. There is only one term in 
the sum for the Pl/2 wave. 

For the D waves a simple polynomial is used 

S2± ^ qHc± + d±q^ + e^q") (5) 

and only the lowest order was used for the F waves 

^3± = //± (6) 
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The r]e± were parameterized allowing a different threshold, qr, for each partial wave. Below the 
threshold the values were unity and above threshold the form 



(7) 



was used except for the D3/2 partial wave which was represented by 



V ■ 



.i=l 



_q_ 



(8) 



The highest power of j varied with partial wave, the greatest being 3 in the case of the P3/2 wave. 
For the F7/2 wave 77 was taken as unity. The values of the parameters for our best fit are given in 
Tables U and ini 

The normalizations of 55 of the data sets was varied in the fitting procedure and their distribution 
is given in Fig. [T] The standard deviation of the normalizations was 2.2%, less than the 3% used for 
AN in Eq. [TJ The mean normalization was 1.0004 and only two normalizations gave an adjustment 
greater than 5%. The total is 2051 for a per data point of 1.09. 
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FIG. 2: T=l phase shifts and 77's obtained in the present work (solid line) compared with VPI (dotted 
line), Hashimoto (dots) and Martin (dashed line). The vertical lines show the thresholds; dash-dot: pion 
production threshold, long-dashed: threshold for production of a theta particle and a pion, dotted: threshold 
for production of a delta and a AT"'', and short-dashed: production by a K*(892) nucleon. 



The phase shifts for the best fit are shown in Fig. [2] compared with three other anal yses discussed 
below. At low energy (below threshold) the results are very similar to those of Martin |l3| . 

Aside from the main fit, we performed tests for the degree of uniqueness of the minimum. Values 
of the parameters were altered in varying degrees such that the values were very large (generally 



6 










10 




20 


— 

OJ — 


30 








40 








50 




60 




70 









-5 


of" 


-10 


3 




Kf 


-15 




-20 




-25 
15 








P3/2 


/ ^ 

/ ^ 




D3/£ 























D5/2 





F5/2 











F7/2 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

P,,H (GeV/c) 



0.4 0.6 0.8 1.0 1.2 1.4 1.6 

P,,, (GeV/c) 



1.00 

0.95 

0.90 

0.85 
1.05 

1.00 

■ 0.95 

0.90 



0.85 
1.05 



1.00 
■ 0.95 
0.90 




D5/2 




F5/2 



r F7/2 ; 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 "'^0.4 0.6 0.8 1.0 1.2 1.4 1.6 
P,.„. (GeV/c) P,,, (GeV/c) 



FIG. 3: T=l phase shifts and ?7's obtained in the present work (sohd line) compared with the variations. The 
sohd hne is the principal fit, the dashed lines correspond to fits a, b, c, e, and f, and the dash-dot line to the 
poorer fit d. The designations of the cases correspond to those in Table HITl 



of the order of 10,000) and the minimum search was restarted. This was done 6 times and the values 
of at the minima found are given in table IIIII The values of were about 9 above the best fit 
except for case d where was about 20 above the best fit. 

The fits found at these other minima are very similar to the original fit. The comparison is shown in 
Fig. [3l There it is seen that the fits are nearly identical below threshold (with the possible exception 
of case d) but above threshold there is a noticeable variation. In principle one can eliminate all 
fits but the principal solution on the grounds that the difference in is considerably larger than 
unity. In practice, however, since one can never be sure that the "true" minimum has been found, we 
consider the variation among the fits (excluding fit d) as a conservative estimate of the error in the 
determination of the phase shifts and inelasticities. 



Other Work for T=l 



Leaving aside many of the earlier analyses [16|, [171, [181, [23, |2J, [251, [261, [271, [201 there are three modern 
analyses with which we compare. 



1. Analysis of Hashimoto 

Hashimoto [l5l| performed an energy-independent analysis for momenta from 0.6 to 1.5 GeV/c. Such 
an analysis has the advantage that no theoretical prejudices are inserted in the parameterization of the 
energy dependence but the disadvantage that S-matrix elements obtained at one energy do not share 
information from nearby energies which can result in large fluctuations from one energy to another. 
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FIG. 4: Top: Contribution of the various partial waves to the reaction cross section. Bottom: Comparison of 
the reaction cross section from the present work with that of the VPI group and Martin. 
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Case 




xVn 


Scatt. Len. (fm) 


Pl/2 (fm^) 


P3/2 (fm^) 


Best Fit 


2031.05 


1.080 


-0.308 


-0.092 


0.103 


a 


2042.47 


1.086 


-0.311 


-0.055 


0.046 


b 


2039.94 


1.085 


-0.310 


-0.067 


0.058 


c 


2040.00 


1.085 


-0.310 


-0.066 


0.057 


d 


2051.98 


1.091 


-0.313 


-0.029 


0.009 


e 


2039.67 


1.085 


-0.308 


-0.084 


0.094 


f 


2039.02 


1.085 


-0.311 


-0.060 


0.051 



TABLE III: Values of for the best fit and the variations made in the present work for T=l. The columns 
labeled Pl/2 and P3/2 contain scattering volumes. 

Observables sometimes have to be grouped to have enough data at a given energy. When amphtudes 
are slowly changing this does not lead to problems but at the threshold for pion production, where 
some of the amplitudes change rapidly, it can. Structure was seen in this analysis in the P3/2 and 
other partial waves. 

2. Analysis of the VPI Group 

The VPI group [l^, [2l|, [2^ has published three analyses. In these fits the S-matrix elements for 
each partial wave were parameterized in the form of a K-matrix with one elastic and one inelastic 
channel. Structure is seen in these fits in the P3/2 wave very similar to that of Hashimoto. This led 
the group to suggest that there was a resonance in this wave and two analyses give the mass pole at 
1.780 and 1.796 GeV/c^, corresponding to beam momenta of 0.971 and 1.005 GeV/c. 

The onset of one-pion inelasticity was assumed to come about by the production of intermediate 
particles which then decayed to the KNtt system. Such a model is made very plausible by the 
experimental fact that pion production docs not start at its threshold but some 200 MeV/c higher in 
momentum (see Figs. |3]and[5|). Since the A threshold is 340 MeV/c above pion threshold and the 
delta has considerable width, so that the production can start below that, it was the prime candidate. 
The if* (892) was considered as well but it has a higher threshold. 

The VPI work [2l|, [g^l allowed renormalization in much the same way as described above in the 
present work. However, unlike the present case, substantial renormalization occurred with a number 
of the data sets being renormalized by more than 5%. 

We inserted the VPI solution (including the G and H waves) into our program to compare with 
the data base used here to find a of 2810 for the 1880 points (a ratio of 1.49), which implies a fit 
very similar in quality to the original fit which was 1.73 or 1.25 depending on the case. Two points 
in the threshold of the reaction cross section (at 0.735 and 0.785 GeV/c) contribute 270 to the x^- If 
we remove these two points the ratio drops to 1.35. One cannot expect a closer agreement since their 
parameters were not fit to our data base and the phase shifts taken from the paper [2^ (or the web 
site) are quoted with a limited accuracy. 

3. Analysis of Martin 

Martin performed an analysis in which the partial-wave amplitudes were parameterized directly. 
The result is equivalent to the partial-wave expansion in terms of Se± and rii± discussed above and 
the correspondence is easily made. 
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FIG. 5: Partitioning of the T=l reaction cross sections among partial waves for Hyslop et al. |22j and Martin 
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FIG. 6: Argand plot of the T=l partial waves. 



Suppressing the partial-wave index, the Martin amplitudes have the form 

Ciq) 



F 



A{q)+l[l + 9{q-qo)RHq)]C{q) 



(9) 



^2f + l 



1 + iq/qo) 



2e+i ■ 



(10) 



A{q) was approximated by a polynomial in q/qm and R{q) a polynomial in ^ ^ ^ where qo is the 

reaction (pion production) threshold (0.3106 GeV/c in the center of mass) and q^, is the maximum 
C. M. momentum considered in the fit. Thus, the inelasticity (expressed in terms of R) has a form 
similar to that used in the present work but only a single threshold (qo) was used for all partial waves. 

In his fit Martin used G and H waves which were estimated by means of a dispersion relation 
calculation [s^l • He did not give these values and they are no longer available so we were forced to set 
them to zero. These waves are only important at the higher energies so did not affect the calculations 
done for the nuclear scattering P, However, for momenta above ~ 1.2 GeV/c the Kaon-nucleon 
amplitude can not be accurately obtained from Martin's phase shifts alone. 

Martin listed several re-normalizations chosen by hand. The average value of these re-normalizations 
was 0.98 and the standard deviation was 0.067. Since these normalizations were changed by hand, no 
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FIG. 7; r=l P3/2 phase shift in the region just above threshold. The vertical line marks the A threshold. 
The chain-dash curve, taken from the paper by Wyborny et al. [s^l and is the result of a calculation of the 
final-state interaction between the and the A. The other curves have the same meaning as in Fig. [2] 



was accorded to them by Martin. With a fixed error of 0.03 they would have contributed 64.7 to 
the x^- 

The Martin solution was also compared with our data base and we found a of 3023 for a ratio 
of 1.61. The same two reaction cross section points which contributed large in the case of the 
VPI group give a contribution to the of 446. If we remove these two points, the ratio becomes 
1.37. Since the normalizations for the data were refit the test is not completely valid. The major 
discrepancy occurs for large beam momenta where the lack of G and H partial waves used in the 
original fit is most important. Removing the two reaction cross section points and all data above 1.2 
GeV/c, we find a x^ of 1372 for the 1132 points of the reduced data set or a ratio of 1.21 so that the 
quality of the fit approaches that quoted in the original paper which was 1.08. 



IV. DISCUSSION FOR T=l 



The values of the phase shifts and inelasticities for the obtained partial waves are given, compared 
with other analyses, in Fig. [21 The behavior of the amplitudes is also given in an Argand diagram in 
Fig. [SI We note that there is no counter clockwise behavior in the P3/2 partial wave as was observed 
in Nakajima et al. 2b\ or Arndt et al. [20l [2l|. 



A. Reaction Cross Section 



The partitioning of the reaction cross section among the partial waves is given in Fig. [H The lower 
part of Fig. [3| gives the comparison of the reaction cross sections calculated with those from Arndt et 
al. [10] and Martin [l^. One notices that just above the threshold for pion production the reaction 
cross section remains very small for about 200 MeV/c. A similar thing is seen in nucleon-nucleon 
scattering and the usual interpretation is that pion production proceeds primarily by A production 
with the threshold for the reaction cross section being governed by the mass of the A. 

In the present reaction, the threshold behavior can perhaps be explained again by A production 
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FIG. 8: Diagrams for pion production. Compound objects are indicated by the double line. 

although the exchange of a single pion to form the A is not possible since a pion cannot couple to 
the kaon. Another possibility is for the reaction to proceed by the excitation of the if* (892). The 
VPI group ^] used these two intermediate states to model the threshold behavior. As can be seen 
from the insert in Fig. U] they obtain a rise in the reaction cross section very similar to that obtained 
by Martin, neither of which is sufficiently rapid. In the case of the VPI group [21] this is because 
the Jr*(892) has a higher threshold (1.077 GeV/c in the laboratory) and, while the A has about the 
right threshold (0.886 GeV/c in the Lab), it has a significant width which causes the onset of the 
reaction cross section to be gradual. However, Oset and Vicente Vacas [5^ had a reasonable success 
in reproducing the data by combining this mechanism with a chiral symmetry calculation. 

Wyborny et al. (56l | included final-state interactions between the A and kaon in the pion production 
channel and were able to produce a maximum in the P3/2 phase shift very similar to the one observed. 
Their result gives a shape closer to our results (or Martin's) than to those of Hashimoto [l3| or Watts 
et al. [28| with which they compared. Figure [7| shows the comparison of their calculation with the 
four sets of phase shifts considered here. 

In order to have a rapid onset for the reaction cross section through the mechanism of an interme- 
diate particle production, that particle must have a small width. An interesting possibility is pion 
production by the intermediate step of K + N ^ +'k where 9^ is the much discussed pentaquark. 
Hyodo et al. [STj ] have considered pion production in this reaction by this mechanism. The mass of the 
proposed pentaquark is thought to be around 1.54 GeV/c^ and its width very small. The threshold 
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for production by this means (0.758 GeV/c in the Lab) is shghtly lower than the A threshold. The 
diagrams for these three possibilities are given in Fig. [8] 

The separation of the reaction cross section into its partial- wave components may be of some help in 
sorting out the reaction mechanisms (see Figs. |3]and[n]). If one assumes that the intermediate particle 
is produced at rest in the center of mass then it should be in a relative s-wave. For the case of the A 
(with spin-parity 3/2-I-) in an s-wave, the only incident partial wave allowed is D3/2. It is interesting 
to note that in the present work, and in Martin's analysis, this partial wave has very little reaction 
cross section and the VPI analysis [2l| has only a moderate contribution. One possible explanation 
is that the production mechanism does not allow the formation of the A-kaon final state in a relative 
s-wave. If we assume that it is produced in a relative p-wave then possibilities for the initial partial 
wave are -Pl/2, P3/2 and F5/2. We see that the Pl/2 partial wave contributes a significant reaction 
cross section and the i-'3/2 partial wave also contains strength, although at slightly higher energies. 

It is interesting to note that the _D5/2 partial wave is important and, at least in the present work, 
dominates in the region 0.8 to 1.0 GeV/c. The Martin analysis also gives the dominant waves to be 
^3/2, Pl/2 and £'5/2 (Fig. [T5|) . In our work it seems to be the D5/2 partial wave which accounts for 
the rapid rise in the region 0.8 to 0.9 GeV/c. If one were to assume that this increase corresponds to 
the threshold for 0+ production in the s-wave then the 6^ would have to have spin-parity 5/2+. 

Bland et al. [13] measured the angular dependence of the pion production directly and found that A 
production is in the Pl/2 and P3/2 waves. They attempted a partial- wave analysis of the production 
although it was only possible to include a limited number of amplitudes. Their analysis indicated that 
the two P waves contributed about equally. They did not include the _D5/2 wave in their study. In 
this same work ^\ a model of A formation by p exchange was also presented which indicated that 
the expected ratio of P3/2 to Pl/2 formation was 5:1, similar to what we find (Fig. [3]). 

Bland et al. (47j also found strong evidence for interference with an even parity partial wave. This 
effect became much larger for low beam momenta and invariant masses away from the central mass 
of the A. This even partial wave was only present in the final charge states K°p7r+ and K+n7r+ 
and not in the state iT^pTr" where the kaon-nucleon system has isospin unity. Also observed in this 
paper was a strong asymmetry in the Daliz plot which, again, was not seen in the JT+pTr^ final state. 
These observations suggest that the asymmetry may be linked to the production of a T=0 KN final 
state, again consistent with the presence of an isospin particle. They [13] also studied the K*{%%2) 
production and were able to say that it occurred in a low angular momentum state consistent with 
our observation of a rise in the s-wave contribution in the region of the K*{%92) threshold. 

B. Scattering Length 

The value of the T=l scattering length has been very stable over the years. Levy-Leblond and 
Gourd in |58| obtained a scattering length of -0.34 fm in what was probably the first analysis. Hyslop 
et al. [22| obtained a value of -0.33 fm, Cutkosky et al. [27] found ~0.28± 0.06 fm and Martin \W{ 
found -0.32 fm. 

We obtain a central value of the s-wave scattering length of -0.308 fm from the best fit. One can 
estimate from the variation of minima in Table Hill an error of ± 0.002 fm. Fixing the scattering length 
at various values and re-fitting the rest of the parameters, one finds again an error of ± 0.002 fm. 
Adding the two errors in quadrature we can quote a value of -0.308 ± 0.003 fm. While this error is 
that obtained with the present data set, it is so small that one would have to expect that the addition 
of new data would lead to a change of the same order or larger. We note that the average from Refs. 
[2^ . [27j and [l^ is -0.031 fm so our central value is what might be expected. 

On the theory side, Barnes and Swanson ^ obtained, in a quark Born approximation, a scattering 
length of -0.35 fm. When this value was corrected for unitarity, by solving for scattering from a 
potential, the value became -0.22 fm. 
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C. Scattering Volumes 

The scattering volumes have always been estimated to be very small and are poorly determined. In 
early work, the p-waves were often neglected as input to the T=0 determinations. Cutkosky et al. [2^ 
obtained -0.04± 0.03 fm^ for the Pl/2 wave and 0.02 ± 0.02 fm^ for the P3/2 wave, for example. The 
values obtained in the present work are very dependent on the minimum found, unlike the scattering 
len gth. For the higher values of the scattering volumes are small, in agreement with Cutkosky et 
al. [27| but for the best fit, the volumes are somewhat larger in magnitude. 
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FIG. 9: Distribution of the values of for isospin for all data. 



V. ANALYSIS FOR r=0 
A. Data Treatment 

The treatment of the isospin zero amplitude is more difficult since there exist no free neutron 
targets. What is normally used is scattering from the neutron contained in the deuteron. Modern 
analyses 'so', '60'| suggest that scattering from the meson cloud is not a problem. There are only 
slightly more than half of the number of data points compared to the T—1 case and the reliability 
of the data is less, given that it must be extracted from deuteron data with corrections. Again most 
of the data can be obtained from compilations fso', 'Sl', but we try to give references to original 
papers [28', i4l, [H, HI [Ml, IsS, il, [64, 65, 66, 61, 68, 69, 70, [Zll, III- Normally the neutron result is 
obtained assuming that the interactions of the kaon with the proton and neutron are independent. 
This is usually a good approximation except at low energies. However, the T=l reaction dominates 
the deuteron (for example the total cross section on the deuteron is |cr(T = 1) + \(y{T = 0)) so a 
good knowledge of the T—\ amplitudes is needed. We used T=l amplitudes from out best fit as input 
to the T=0 determination. 

The data used in our standard fit are all from incoherent deuteron scattering and deuteron 
total cross sections. There are 92 points of total cross section data, 17 points of one and two pion 
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FIG. 10: Partitioning of the r=0 reaction cross sections among partial waves for the present work with the 
standard reaction data. 



production data, 336 points of elastic cross sections from the neutron, 43 points of elastic polarization 
from the neutron, 657 points of charge exchange cross section data and 44 points of charge exchange 
polarization data for a total of 1189 points. 

As might be expected, the pruning of the T—Q data is more significant and problematic than the 
r=l case. In order to choose which points to eliminate, a preliminary fit was made with all data 
included. By binning the values of (bins of unit size were chosen) a distribution was obtained. 
This distribution, along with the expected distribution obtained from a distribution with one degree 
of freedom is shown in Fig. [51 We see that the expected number of counts in a bin of unit size falls 
to 0.1 around x^=15. For this reason, the 19 points with > 15 were removed leaving 1170 data 
points to be fit. It is clear that there is still an excess of points below — 15 but it is impossible to 
tell which points to remove so no further pruning was done and we must expect to have a higher x^- 
However, we did make one run with all points giving a x^ > 10 being removed. The results of this fit 
are shown in Table IIVI 



B. Double Scattering 

At very low energies, double scattering in the deuteron gives an important contribution to the total 
cross section. Other observables are not sufficiently well measured that it will have a significant effect. 
The double scattering amplitude at 0° is given by [ill 

/D(e-0)-/z5(k,k) = — y — ^— — z(|k-q|) (11) 



where k is the initial and final (on-shell) momentum of the scattering meson, 2;(p) is the two-body 
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FIG. 11: Partitioning of the T=0 reaction cross sections among partial waves for the present work with 
alternate reaction data. 



form factor, and /(k, q) and /(q, k) are half-off-shell basic scattering amplitudes. 
For the scattering amplitudes we write the off-shell dependence as 

/(q,q') = fov{Q)v{q'); v{k) = l where the form v{q) = ( ^2 ^ V ^^'^^ 



will be assumed. 



M. = 0) = |/d^^,/.(q,k)/„(k,q>(|k-q|) + ^P / ^^.^-^^/^(q.^^^C^, C,) ^ _ ^^3^ 

We will consider only the s- and p-waves for this correction. The isospin 1 s-wave is the strongest 
so we consider the double scattering between it and the s- and p-waves of the neutron. With these 
assumptions 
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(14) 



where x is the cosine of the angle between q and k. The double scattering contribution to the total 
cross section will be 
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FIG. 12; Distribution of the normalization factors for isospin 0. 



^fUl / dVL^x^ (|k - q|: 



nk 



ik-qi: 



(15) 



Where the factor of 2 comes from the fact that there are two orders of scattering possible, 
include charge-exchange scattering we can replace 



To 



fP fn 
JQ Jo 



(16) 



Near zero energy the amplitudes become real (so the principal value terms become very small) and 
the p-wave amplitude goes to zero so that only the first term remains and gives a contribution of 
Stt/q/q. Since the fourth term is small for both reasons we have neglected it. The function z was 
computed using the one-pion-exchange deuteron wave function J 4] which gives a good representation 
of the momentum distribution in the deuteron [75| . While we believe that this correction is needed 
to get a proper fit, the results with it being left out are shown in Table [TVl 



C. Reaction Cross Sections 

Aside from the values given by Hirata et al. [U , the principal reaction cross sections (one and two 
pion production) are from Giacomelli et al. (6^ . They used their deuteron and proton pion production 
data to extract the neutron data. A more consistent way for us is perhaps to use our fit to the proton 
results with their deuteron production data to obtain directly the T=0 reaction cross section. These 
can be obtained directly in the single scattering impulse approximation from the following equations. 



CT/j(T = 1) = crfl(proton); (7^; (deuteron) = crfl(proton) 



(Tfl (neutron) = 2'^r{^ 



1) + 2'^fi(T = 0) 
(17) 
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Case 


v2 

A 


TV 


A. / ' 


Scatt. Len. (fm) 


Pl/2 (fm^) 


P3/2 (fm^) 


Best 


1670.561 


1170 


1.428 


-0.1048 


0.183 


-0.029 


a 


1671.214 


1170 


1.428 


-0.1015 


0.183 


-0.028 


b 


1670.561 


1170 


1.428 


-0.1053 


0.182 


-0.029 


c 


1670.758 


1170 


1.428 


-0.1055 


0.182 


-0.029 


d 


1671.085 


1170 


1.428 


-0.1055 


0.183 


-0.030 


e 


1675.274 


1170 


1.432 


-0.1067 


0.181 


-0.022 


f 


1675.301 


1170 


1.432 


-0.1003 


0.181 


-0.022 


x' < 10 


1540.480 


1162 


1.326 


-0.1027 


0.182 


-0.028 


No Krauss/ Weiss 


1649.073 


1155 


1.428 


-0.0988 


0.177 


-0.027 


No Double 


1662.132 


1170 


1.421 


-0.1166 


0.200 


-0.038 


Alter. Reaction 


1665.105 


1170 


1.423 


-0.0957 


0.183 


-0.021 


Stenger Full 


1714.224 


1190 


1.441 


-0.1096 


0.174 


-0.036 


Stenger Partial 


1692.632 


1186 


1.427 


-0.1036 


0.181 


-0.029 


No Damerell 


1284.174 


1008 


1.274 


-0.0997 


0.205 


-0.025 


KOp 


1895.322 


1319 


1.437 


-0.1069 


0.173 


-0.032 



TABLE IV: Values of for the best fit and the variations made in the present work for r=0. Cases a-f show 
different fits obtained with the full code. Case "No Krauss/ Weiss" corresponds to a fit in which the data of 
Refs. [3 and d were left out. In case "No Double" the double scattering correction was omitted and in case 
"Alter. Reaction" the reaction data was replaced by the average as explained in the text. The case "Stenger 
Full" gives the results with the full Stenger data [76ij included and case "Stenger Partial" the results where 
the most forward points were omitted. For the case "No Damerell" the Damerell et al. data [t^I were omitted 
from the fit. For the line labeled "KOp" the data of Armitage et al. [ttI were included in the fit. 



LJ 


qa {GeV/c) 


71 


72 


73 


Equation 


Sl/2 


0.4350 


3.581 


-4.349 




El 


Pl/2 


0.2515 


-0.689 


0.822 




El 


P3/2 


0.5303 


-9.101 


46.593 


-79.735 


El 


L»3/2 


0.3351 


0.887 






El 


L>5/2 


0.7120 








El 


F5/2 


0.6352 


1.485 






m 



TABLE V: Parameters for the representation of the inelasticity, rjt^j for T=0 using the form of Eq. El for all 
the partial waves except the -F5/2 wave which uses Eq. [71 The threshold for the 1)5/2 wave is Pi,a()=l-55 
GeV/c so the rj for this partial wave can be taken as unity over our fitted range. . 

With not much to choose between the two one might hope to improve the errors by averaging the 
two determinations. Figures [TUl and [TTI show the result of a fit with the original reaction data and the 
alternate reaction data. Table IIVI shows the effect of using this alternative reaction data. 

D. Normalizations 

We maintained the constant error of 3% for most of the normalizations although, for the reasons 
mentioned above, we might expect it to be larger. Figure [T^ shows the distribution of normalizations. 
The distribution is wider than for the T=l case and there is a slight bias. The calculated width of 
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LJ 


a (GeV/c)-(^^+i) 


(GeV/c)-2 


b2,tt (GeV/c)-* 




Equation 


51/2 


-0.531 


-1.206 


1.362 




HI 


Pl/2 


23.765 


3.690 






HI 


P3/2 


-3.808 


2.919 


-10.042 


212.021 


HI 


LJ 


c± (GeV/c)"'' 


d± (GeV/c)~^ 


e± (GeV/c)-'' 




Equation 


D3/2 


12.548 


-22.412 






El 


D5/2 


-7.528 


32.262 


-37.422 




m 


LJ 


/± (GeV/c)-^ 








Equation 


F5/2 


2.836 








El 


F7/2 


-1.731 








m 



TABLE VI: Parameters for the representation of the phase shifts for r=0 using the form of Eq. HI for the S 
and P waves, Eq. [5] for the D waves and Eq. [6] for the F waves. 

the distribution is influenced by one outlier from the Damerefl et al. elastic data [72] . 

E. Other Work 

L Analysis of Hashimoto 

Hashimoto [lal observed a resonance-like structure in the T—0 DZ/2 wave which had been seen 
before [H, [26(] . While we see a loop in the Argand diagram, it is less pronounced than what he 
saw. 

2. Analysis of Martin 

Martin's database was somewhat smaller than the one used here. He used, in addition to the direct 
deuteron data, real parts of the amplitude obtained from a dispersion relation analysis using the 
and K~ total cross sections. The dispersion relations were once subtracted with the subtraction point 
was taken at zero energy which means that the scattering lengths were input. He assumed that the 
isospin scattering length was ± 0.04 fm based on previous analyses. His fit resulted in a scattering 
length of -0.035 fm in agreement with his small input value. 

We inserted Martin's solution in our code to compare with our database. Since he had a smaller 
amount of data and the parameters were not fit to the present data base, one cannot expect a reduced 
very close to what he obtained. We find (for 1170 points) a of 2849 or a ratio of 2.44. for 
the full data set. Restricting the comparison to data below 1.2 GeV/c (because of the missing high 
partial-waves problem mentioned before) we find a of 1430 for 827 points for a ratio of 1.73. His 
original fit (number 2) found a of 924.5 for 760 points or 1.22. 

3. Analysis of the VPI Group 

Hyslop et al. use A production as a model for inelasticity in spite of the fact that the K+A 
final state is forbidden in the T=0 channel. Comparing their amplitudes with our data base we find 
a x^ of 2776 on 1170 points for a ratio of 2.37. They found 3181 for 1746 data points for a ratio of 
1.82. Again, we remark that their parameters were not adjusted to our data base and phase shifts of 
only limited accuracy are available so the same value cannot be expected. 
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FIG. 13: Partitioning of the T=0 reaction cross sections among partial waves for Hyslop et al. [52] and Martin 
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FIG. 14: T—Q phase shifts and rj^s obtained in the present work (solid line) compared with VPI 23] (dotted 
line), Hashimoto [Ts!] (dots) and Martin ^19] (dashed line). The vertical lines show the relevant thresholds, 
dash-dot: pion production threshold, dashed: threshold for production of a . 

VI. RESULTS FOR T=0 

Figure [Ml shows the phase shifts obtained for the T=Q fit compared with three other analyses and 
Tables fVl and IVTl show the parameters which give the phase shifts and inelasticities. Figure [T51 shows 
the behavior of the T=0 phase shifts in Argand plots. There is some structure in the plots but nothing 
that can surely be associated with resonances. 

A. Reaction Channels 

The isospin channel is quite interesting from the point of view of the mechanism for pion pro- 
duction. Because of isospin, A production is not allowed, nor is the production of the isoscalar 
Indeed, the reaction cross section is seen to be smooth in the region of these two thresholds, unlike 
the T=l case. The K*(892) production is permitted and one does observe a rapid rise around where 
it would be expected (Figs. [TOl and [TT|) . This is particularly true of our results but is consistent with 
the reaction cross sections of the VPI group and Martin as is seen in Fig. [T31 Since the K*{892) has 
spin-parity 1^, if it is produced with the nucleon in a relative s-state, the spin-parity values possible 
are l/2~ and 3/2~ which correspond to the incident waves of 51/2 and D3/2. Indeed, we see the 51/2 
wave giving an important contribution to the reaction cross section at the threshold and the D3/2 
wave is the largest single contributor. Somewhat surprising is the dominance (or at least importance) 
of the D3/2 channel. This is the partial wave expected for the (isospin forbidden) production of a A 
in the s-wave. However, this wave may be simply the dominant wave for the non-isobar production 
as well as receiving a contribution for the if* (892) production. 
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FIG. 15: 



Argand plots of the real and imaginary parts of for r=0. 



B. Scattering Length 



The r=0 scattering length has a rather checkered history. Levy-Leblond and Gourdin [58| obtained 
a value of -0.05 fm, with large but unspecified errors. Stenger et al. [t^ found a value of +0.04 
fm (a value which was commonly used in dispersion relation work There were also suggestions 

that it might be positive and large Q • Presumably based on this previous work, Martin [l^] set the 
scattering length to zero with an error of 0.04 fm as the subtraction point for his dispersion relation 
constraint. Martin's fit gives -0.035 fm, although he states that it has a large error. Later Martin 
gives [3| a value of 0.02 fm and then in still later work [2^ he found -0.23±0.18 fm. 

Barnes and Swanson [9] obtained a theoretical scattering length of -0.12 fm from a Born quark 
model. To compare with the experimental value they performed their own extrapolation to zero 
energy based on single-energy analyses and found -0.09 or -0.17 fm depending on which analysis they 
used. 

From Hyslop et al. we extrapolate a value of -0.019 although they quote in the paper a value of 
zero. It can be seen from Fig. [T4]that the behavior of the isospin s-wave in their fit is rather different 
from ours and has a great deal of variation in the low-energy region where there are no data. It can 
only be assumed that the variation in this case is a result of a fit at higher energies. It is seen that 
the trend of the curve above 1.1 GeV/c is noticeably different from the other determinations. 
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FIG. 16: T=0 phase shifts and ?7's obtained in the present work (solid line) compared with the variations. 
The sohd line is the best fit, the dashed lines correspond to fits a, b, c, and d, and the dash-dot lines to the 
poorer fits e and f where the labeling of the cases are indicated in Table IIVI 

For the T=Q s-wave scattering length we adopt for the central value our best fit value of -0.105 fm 
from Table HVl The error for the uncertainty in the minimum is estimated from Table ITVl to be 0.002 
fm. In order to estimate the statistical error the fit was redone for several fixed values of the scattering 
length varying all other parameters. From the resulting curve the error can be estimated to be 
0.01 fm. For values of scattering length close to the central value, the curve is symmetric but for 
larg er deviations it is not, rising steeply for small values. For the value of -0.035 fm given by Martin 
[l9|, the corresponds to 8.6 standard deviations from our central value. Since the statistical error 
dominates, we take the scattering length as —0.105 ± 0.01 fm. This error does not include possible 
errors from the variations in the database. Observations on the change in value from the omission of 
data sets as shown in Table IIVI suggest that the error from this source can be expected to be of the 
same order or slightly smaller. 

C. Scattering Volumes 

The Pl/2 scattering volume is well determined by the fit to be 0.183 ± 0.005 fm where the error 
comes from an examination of Table IIVI assuming that there is no reason to exclude the double 
scattering correction or the Damerell data. 

The P3/2 scattering volume is smaller and more poorly determined. Again from Table ITVl we take 
the value of -0.029 ± 0.008 fm. 
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FIG. 17: Square-well spin-orbit model. The solid line gives the results of the present amplitude analysis, the 
dotted lines the model described in the text without a central term and the dash-dot lines the model with the 
central potential which gives the correct s-wave scattering length. 



D. Spin-orbit Splitting 

The phase shifts for ^ > shown in Fig. [T3] display a remarkable symmetry below threshold. All 
phase shifts for j = i — ^ are positive and all of those for j = £ + are negative. In order to see 
how far a pure spin-orbit interaction would go toward explaining this behavior, we calculated a simple 
model consisting of scattering from a square well potential with strength V where 

The radius of the well was taken to be i? = 0.85 fm (corresponding to an rms radius of 0.66 fm). and 
the strength, Vq was chosen to be 0.36 GeV. The results are shown in Fig. [T7]with the dashed hue. 
The angular momentum barrier changes a great deal from one value of ^ to another and the potential 
strength also changes over a significant range. The rather remarkable agreement indicates that for 
£ > the phase shifts in the lower energy region are described by a pure spin-orbit interaction. 

Such a potential gives no contribution to the s-wave. If we introduce a central potential (independent 
of t) in all partial waves of strength 0.04 GeV we obtain an s-wave scattering length of -0.11 fm (in 
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agreement with our determination from the data). The result for the higher partial waves is shown 
in Fig. [17] by the dash-dot lines. Thus, including a central potential of sufhcient strength to give the 
moderate s-wave scattering length does not destroy the good agreement seen before. 



E. Variations in the fit 



As in the T=l case, a number of different minima were found corresponding to different starting 
points. The basic properties of the different fits are given in Table ITVl (caes a-f) and the variations of 
the phase shifts are shown in Fig. [16 

The elastic data by Stenger et al. 76| was not included in the general fit. These data were among 
the first to find a very small scattering length for the T—0 channel. They consist of charge exchange 
and elastic scattering from the deuteron. In the elastic scattering there was no separation of scattering 
from the neutron or proton or, indeed, coherently from the deuteron (leaving it intact). Thus, for 
the elastic scattering the cross section from the proton should be added to that of the neutron and 
coherence must be taken into account as well. The contamination from coherent scattering is largest 
for small angles. For the charge-exchange cross section small angle scattering (small momentum 
transfer) tends (without spin-flip) to leave two protons in a triplet s-wave state. Since this state is 
blocked by the Pauli principle and the spin-flip amplitude is small at small angles, this effect leads to 
a very large suppression of the charge- exchange cross section such that it is far from charge exchange 
on a free neutron. The data were taken in large angle bins (0.4 in cos 9). In order to estimate the 
effect of leaving out this data set we made two runs, one in which the two lowest energies were fully 
included in the fit and one in which the most forward points (at cos9 = 0.8) were excluded. The 
results are summarized in Table ITVl 

There are two modern data sets of total cross section data, those of Krauss et al. fS] and Weiss et 
al. ^ . In the fitting process these data suffer a renormalization (down) of about 4% which is greater 
than might be expected. To see the effect of these data on the fit, a run was made with them left out. 
The result is shown in Table IIVI 

Martin [l^ comments that the Damerell et al. [T^I data fit poorly with the rest of the data base. 
For this reason we made a fit with this data removed. Since there is no a priori reason to mistrust 
these data they were used in all of the other fits. The results of this fit are shown in Table IIVI 
The normalization of the elastic scattering (K+n) is greater than 20% and there is some shift in the 
low-energy parameters obtained. 

Two experiments have been performed [t^ . [TSj using the inverse charge-exchange reaction on proton 
targets with JT-long beams . It has been recognized for some time that there is great difficulty in 
controlling the normalization of these beams because of regeneration of K-shoTt mesons. For this 
reason we did not use these data in the principal fit. However, we did include one of the sets [T^I in a 
run to see the possible effect. The results are shown in Table ITVl Table rvTIl gives the normalization 
factors which result from the fit. The normalization errors in the fit were taken as 0.1 through 0.95 
GeV/c and 0.2 above that as suggested in the experimental paper fn\ . 



VII. SUMMARY AND CONCLUSIONS 



We have presented an easy-to-use parameterization of the K'^ Nucleon amplitudes. As mentioned 
in the introduction, there is a need for reliable amplitudes for several purposes. 

Our representation for s- and p-waves is equivalent to the effective-range expansion which, for 
s- waves, reads 



kcotS = l/a+ -rok'^ . . . (19) 

This form was derived originally for a potential interaction but was shown to be valid for an effective 
field theory [79l . [80| . More recently it has been shown to arise from renormalization group calculations 
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U. 1 \JiJ 


0.75 


1.025 


0.85 


0.684 


0.95 


0.667 


1.05 


0.695 


1.15 


0.881 


1.25 


1.082 


1.35 


1.210 


1.45 


1.356 



TABLE VII: Normalization factors for the data of Armitage et al. [ttI ]. 

[sH ]. Since one would tend to believe that this is an appropriate expansion it would seem to be ill 
advised to set the scattering length to zero. 

We now discuss several points which are particularly interesting for their physics potential. 

1) The ratio of the scattering length for T=\ to that of T'=0 may be an interesting quantity. For 
example, in the work of Barnes and Swanson Q this ratio depends only on the ratio of the sizes of the 
kaon and nucleon and the strange quark mass, being independent of the absolute size of the hadronic 
systems. While this is only approximately true in their work (i.e. only in Born approximation), it 
indicates that this may be a quantity which is sensitive to only a restricted set of physical parameters. 
Our value for this ratio is 2.9 ± 0.3. 

2) The reaction cross section in the T={) state may provide an interesting piece of data for the cal- 
culation of pion production. Here the usual dominant mechanism (A production) is isospin forbidden 
so that other mechanisms will be more apparent. Calculations of the type of Oset and Vicente Vacas 
[ssj might be interesting for the T=0 channel. 

3) The simple form of the ^ > phase shifts for the T=0 amplitude is remarkable. To date, no 
theoretical model has been able to reproduce this feature althugh the work of Biittgen et al. [H^ was 
able to get a moderately good representation of the data with some degree of phenomenology. 

4) The presence of a narrow pentaquark state would facilitate the understanding of pion production 
in the T=l channel. In fact, the best way to look for a narrow resonance may be to produce it and 
look for a sudden change in the inelastic cross section. This was the way in which the existence of 
the J/^* was first indicated. If this is indeed the explanation of the rapid rise in pion production and 
if the final — tt state has relative angular momentum zero, then the spin-parity of the Q must be 
D5/2+. Thus, the resonance, as seen directly in the T=0 channel, would be in the F5/2 partial wave. 
We have seen that this partial wave is very attractive with a large angular momentum barrier so that 
a narrow "molecular" state is possible. 

We thank Jean-Pierre Dedonder for comments after a careful reading of the manuscript. 
This work was supported by the National Science Foundation under contract PIIY-0099729. 
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